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Abstract 

We present and analyze a simplified stochastic aJ7— dynamo model which is designed to 
assess the influence of additive and multiplicative noises, non-normality of dynamo equation, and 
nonlinearity of the a— effect and turbulent diffusivity, on the generation of a large-scale magnetic 
field in the subcritical case. Our model incorporates random fluctuations in the a— parameter 
and additive noise arising from the small-scale fluctuations of magnetic and turbulent velocity 
fields. We show that the noise effects along with non-normality can lead to the stochastic 
amplification of the magnetic field even in the subcritical case. The criteria for the stochastic 
instability during the early kinematic stage are established and the critical value for the intensity 
of multiplicative noise due to a— fluctuations is found. We obtain numerical solutions of non- 
linear stochastic differential equations and find the series of phase transitions induced by random 
fluctuations in the a— parameter. 
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1 Introduction 

The understanding of the generation and maintenance of a large scale magnetic field in astro- 
physical objects is a problem of exceptional importance and difficulty. It is widely accepted that 
the magnetic field is generated by the turbulent flow of the electrically-conducting fluid. Inho- 
mogeneous velocity fluctuations stretch magnetic lines and amplify the magnetic fleld. These 
small scale fluctuations of turbulent flow are primarily responsible for the generation of mag- 
netic flelds. The problem is that it is difficult to resolve them. Thus their influence on the 
resolved large-scale magnetic fleld has to be modelled. A traditional closure scheme is based on 
the a— effect according to which small-scale fluctuations can be described by an average term 
involving the curl of the mean magnetic fleld, B, written as V x (aB). The asymptotic analysis 
of the induction equation exploiting the assumption of two separated scales for turbulent flow 
leads to the effective macroscopic equation for the large scale magnetic field B(t,x) [1-5] 



where u is the mean velocity fleld of the turbulent flow, a is the coefficient of the a-effect and 
/3 is the turbulent magnetic diffusivity. Traditionally the phenomenon of generation of magnetic 
fields was analyzed by considering perturbations of the trivial state B = and looking for 
exponential solutions to the deterministic PDE with appropriate boundary conditions. While 
this standard stability analysis successfully predicts the dynamo action for the supercritical case, 
there are situations for which this eigenvalue analysis fails to predict the subcritical onset of 
instability [6-10]. It was pointed out in [9^ that the closure scheme involving only deterministic 
a/3— parameterization is not completely satisfactory since unresolved fluctuations may produce 
random terms on the right hand side of the dynamo equation (^. Moreover, it follows from 
astronomical observations that large scale magnetic fields exhibit a rich random variability both 
in space and time that cannot be described by the deterministic equation 

The importance of noise effects in the dynamo problem has been recognized previously and 
several attempts have been made to account for the effects of spatial and temporal fluctuations 
in small scale magnetic and velocity flelds on the generation process. Kraichnan considered 
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fluctuations in the a-parameter and found a negative contribution to turbulent diffusivity 
from helicity fluctuations 11 . Hoyng with colleagues in [12-14] studied in detail the effect of 
random fluctuations in the a-parameter by considering the system of stochastic linear equations 
for eigenmodes corresponding to the dynamo equations. They found the excitation of those 
modes such that their magnetic energy is proportional to 7^^, where 7 is the damping rate. 
Stochastic dynamics of magnetic field generation have been also analyzed by Farrell and loannou 
in where they examine a mechanism by which small-scale fluctuations excite the large scale 
magnetic field. They modelled these fiuctuations by an additive noise term in the mean field 
equation and identified the crucial role of non-normality on the dynamo process. Numerical 
simulations of the magnetoconvection equations were performed in with analysis of effects 
of noise and non-normal transient growth. Inhomogeneous turbulent helicity fluctuations were 
considered in IF. One should mention the dynamo model that exhibits aperiodic switching 
between regular behavior and chaotic oscillation . Stochastic dynamo theory, using the term 
of an incoherent dynamo, was proposed by Vishniac and Brandenburg in JH]. They showed how 
random fluctuations in the helicity can generate a large-scale magnetic field for the aJl— dynamo. 
We note that this model is closely related to the present work. However, they did not consider 
the transient amplification due to the non-normality of the dynamo equation. It turns out that 
non-normal dynamical systems exhibit an extraordinary sensitivity to random perturbations 
which leads to great amplification of the second moments of the stochastic dynamical systems 
(see j7| 1191 \'2U\). Our recent work has demonstrated the possibility for stochastic magnetic 
energy amplification in the subcritical situation where the dynamo number is less than critical 
jlUj . These observations motivate further studies of noise and nonlinearity effects which we 
consider below in the context of a simplified stochastic no-z model (see j21|). An issue we 
address in this paper is how the random fluctuations may be appropriately incorporated into 
the classical aO— dynamo model. Our analysis especially focuses upon the multiplicative noise 
due to random fiuctuations in the a-parameter and the non-normality of the dynamo equation 
operator. Recall that an operator is said to be non-normal if it does not commute with its 
adjoint in the corresponding scalar product. The determination of the effect of the noise, along 
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with non-normality, on the amphfication of magnetic energy during the early kinematic stage is 
the primary contribution of this paper. 

In section II, we discuss the deterministic aJl— dynamo model, its equilibrium points and 
transient growth effects. In section III, we consider a linear stochastic model for the subcritical 
case. We derive equations for the second moments and find their stationary values. We demon- 
strate the important differences between the non-normal system and the normal one under the 
influence of additive noise. We explore the influence of multiplicative noise due to fluctuations 
in the o-parameter on the amplification of magnetic energy in the subcritical case during the 
kinematic stage. We derive the criteria under which the second moments grow exponentially 
with time (kinematic regime). Finally, in section IV, we perform numerical experiments showing 
that there exists a series of noise-induced phase transitions in the aO— dynamo model. 

2 Stochastic ail— dynamo model in a thin-disk approximation 

Following |H1 11] we consider here the thin-disk approximation to the dynamo equation for spi- 
ral galaxies, in which a turbulent disk of conducting fluid of uniform thickness 2h and ra- 
dius R {R ^ h) rotates with angular velocity Q,{r). We restrict ourselves to the case of 
the afi— dynamo for which the differential rotation dominates over the a-effect. The governing 
equations for the components of the axisymmetric magnetic field, {t, r, z) and {t, r, z) 
in the polar cylindrical coordinates (r, ip, z) can be written as, 

dBr d, , , B'^Br (3d fldBr\ 



dt fJz ' ' dz"^ r dr \r dr 

where = rdQ/dr is the measure of differential rotation (usually rdQ/dr < 0), and frit, z,r) 
and fip{t,z,r) are the stochastic terms describing unresolved turbulent fluctuations. The com- 
ponents Br {t, r, z) and B^p {t, r, z) obey vacuum boundary conditions on the thin disc surfaces 
Br^^ (t, r, -h) = 0, Br,^ (t, r,h) = (see i]). 

In what follows we neglect the spatial structure of the magnetic field along the radius r and 
the height of the galaxy adopting the well-known no-z model I^. The aim here is to concentrate 



on the studies of the influence of random fluctuations and non-normahty on the dynamo process. 
We also take into account the non-hnear backreaction [,V^\ l,S4j . The dynamical system for the 
azimuthal, (t) , and radial, Bj. (t) , components of the magnetic field B can be written then 
in the form 

dB^_ ^^/?(|B|) 

^ - 9^Br - -^^B^ . (3) 
This dynamical system may be regarded as an one-mode approximation of @. It should be 
noted that the deterministic critical conditions for the generation of a magnetic field are the 
same for both models © and ah^(3~'^\gi^\ > vr^/lG jlj. Here we introduce the non-linear 
functions a(|B|) and /3(|B|) describing the quenching of the a— effect and turbulent magnetic 
diffusivity. It is well-known that the current theories disagree about how the dynamo coefficients 
a(|B|) and /9(|B|) are suppressed by the mean field B [231 EH EHl EHl EZl • Here we describe the 
dynamo quenching by using the standard forms [3] 

a{\B\) = ao{l + UB^/B,gfy\ /J(|B|) = /3o (^1 + ^^^^^y^^ , (4) 

where ka and k/^ are constants of order one, and B^q is the equipartition strength. We note that 
for the afi— dynamo the azimuthal component B^ (t) is much larger than the radial field Br (t) , 
therefore, B^ ~ B'^. The issue of the strong dependence of a and f3 on the magnetic Reynolds 
number is still controversial and we did not include it in our analysis. 

Regarding random perturbations in the system @ , we adopt a phenomenological mesoscopic 
approach in which the parameters of dynamo equations are assumed to be random functions of 
time j3U[ I31j . The multiplicative noise £,ait) in @ stands for the rapid random fluctuations in 
the parameter a. One can show that they are more important than the random fluctuations in 
the turbulent magnetic diffusivity /? [I'll 1131 114j . The additive noise term (t) represents the 
stochastic forcing arising from the small-scale fluctuations of magnetic and turbulent velocity 
fields (S]. Both random terms are assumed to be independent Gaussian white noises with zero 
means < ^a(t) >= 0, < (,r(t) >= and correlations: 

<Ut)Ca{s)>=2D^6{t-s), <Ut)Us)>=2Dr5it-s), (5) 



where and are the intensities of the noises and the angular brackets , < • >, denote 
the statistical average. One can show that the additive noise in the second equation in Q 
is less important and can be ignored. In particular, the main contribution to the second mo- 
ment of comes from the additive noise in the first equation (see, for example, [SOD- Of 
course, this paper addresses the over-simplified case of magnetic field generation in galaxies. 
Nonetheless, we present this work as an illustration of the influence the random fluctuations and 
non-normality may play in the dynamo process, and which therefore should be accounted for in 
more complicated dynamo modelling like the stochastic PDE @. 

The dynamical system ((21 is well studied for the case when the noise terms are absent, and 
there has been much progress in the prediction of the growth rates induced by both the a— effect 
and differential rotation and the corresponding speed of magnetic waves !23_. However, there 
is considerably less understanding of generation of a magnetic field in the presence of random 
fluctuations. Although the equations © are a theoretical simplification of what really happens 
in galaxies, we strongly believe that it does provide a useful framework for understanding the 
interaction of stochastic perturbation, non-normality and non-linear effects. 

It is convenient to rewrite the governing equations © in a nondimensional form by using 
an equipartition field strength i?eg, a length /i, and a time Jig ^, where VLq is the typical value of 
angular velocity. In terms of the dimensionless parameters 

the stochastic dynamo equations (jSJ can be written in the form of SDE j3L)l I31j 

dBr = -{5ipa{Bp)B^ + e^p{Bp)Br)dt - V2^^a{B^)B^dWi + V2^dW2, (7) 

dB^ = -{gBr + e^p{B^)B^)dt , (8) 

where Wi and W2 are independent standard Wiener processes, and ai and CJ2 are the noise 
intensities 
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Here we introduce the functions ipa{B^) and ipp{B^) describing non-linear quenching 

Now let us discuss the other parameters in (jT)),®, namely, 5 = Ra/Ru and e = tt'^/AR^^. The 
parameter 6 is the characteristic of relative importance of the w— effect and the q— effect. For the 
ail— dynamo the differential rotation dominates over the a-effect, that is, Ra ^ Ruj- Moreover, 
the diffusion time /i^//? is much larger than the time ^, therefore, both parameters 6 and £ are 
small. Their typical values for spiral galaxies are 0.01 — 0.1 {R^^ = 10 — 100, Ra = 0.1 — 1) It 
turns out that for small values of 6 and e , the linear operator (matrix) in (O,® is a highly non- 
normal one, since ~ 1. This can lead to a large transient growth of the azimuthal component 

(t) in the subcritical case. The nonlinear interactions may lead to a further amplification of 
this small disturbance jlOj . The crucial idea behind subcritical transition is that the a— effect 
or the a;— effect might be relatively weak, but the generation and maintenance of the large scale 
magnetic field is still possible. 

3 Deterministic system 

Let us briefly review the dynamics of the system ((T)),© in the absence of noise terms. It takes 
the form 

dB 

^ = -gBr-eMB^)B^- (11) 

The equilibrium points of the system (fTT|) can be found from 5ipa{Bip)B^ + £ipf^{B^)Br = 0, 
gBr + £ip/3{Bip)B^ = 0. The stationary value of the radial component, Br, can be expressed in 
terms of the azimuthal one, B^, 



B^ = _'MB^)B, 
9 



while B^p is a solution of the equation 



p2 

-^}{B^)B^ + 6ipa{B^)B^ = 0. (13) 



Over a wide range of parameters, this equation might possess five solutions including the trivial 
one = . The other stationary points can be found from the equation 

.2 



-(1 + B'^Yil + kaB'^) - (1 + {kp + l)B'^y = 0. 



(14) 



Since there are multiple stable solutions (see Fig.l for e = 0.1, 6 = 0.01, k^ = kp = 1), the 
system ((TT|) can exhibit hysteresis P^ . 



- stable 
unstable 



0.6 0.7 0.8 0.9 1 



1.2 1.3 1.4 



Fig. 1. The dependence of the stationary points of the deterministic system (|11|) 
on the parameter g (subcritical bifurcation). 

Of course, the stochastic system, in which such multiple stable states can coexist (metasta- 
bility), may exhibit random transitions between these states. In general, the transition rate is 
inversely proportional to the corresponding mean first passage time |3n( 131 j. In the next sections 
we discuss the peculiarities of such transitions with non-normal effects. 

The onset of the instability of the trivial equilibrium state B^ = 0, B^^ = can be obtained 
from a standard linear stability analysis. Linearization gives 

dB 



dt 



AB, 



where 





-e -5 




Br 


A = 




, B = 




_-g -£ _ 




B^ 



(15) 



(16) 
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The matrix A has two eigenvalues 

Xi = -e + y/g6, X2 = -e- y/g5. (17) 

The supercritical instability condition (Ai > 0) can be written as g6 > e^. The suhcritical case 
corresponds to g6 < e^. The main purpose of this paper is to study the stochastic amplification 
of the magnetic field B in the suhcritical case. 

Since 5 and e are small parameters and (7 ~ 1, the matrix A is a highly non-normal one. 
Recall that the matrix A is normal, if AA"^ = A-^A, where denotes the Hermitian transpose, 
otherwise it is non-normal. Even in the suhcritical case when both eigenvalues Ai^2 are negative, 
B^p exhibits a large degree of transient growth before the exponential decay. The azimuthal 
component Bip{t), as a solution of the system (|15|1 with the initial conditions 

Bri0) = -2c^, B^{0) = 0, (18) 

has the form 

B^{t) = c(e^^* - e^2*). (19) 

Thus Bip (t) exhibits large transient growth over a timescale of order ^^^^^^^ In ^ before decaying 
exponentially. In other words the transient growth causes a temporary exit from the basin of 
attraction of the linearly stable solution (0,0). In Fig. 2 we plot phase portraits of the system 
jni) for e = 0.1, 5 = 0.01, ka = kp = 1. 




(a) g = 0.8 one stable stationary point (0, 0) (b) g = 0.95 three stable stationary points 
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(subcritical case) 



"-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 

(c) g = 1.1 two stable stationary points (supercritical case) 
Fig. 2. Phase portraits of deterministic system (|11|) 

Similar deterministic low-dimensional models like Hll|) with transient growth have been pro- 
posed in |241 1251 1261 127j to explain the subcritical transition in the Navier-Stokes equations j2H] • 
It will be interesting to analyze the so-called self-killing and self-creating dynamos by Fuchs, 
Radler, and Rheinhard j^H] in terms of non-linear subcritical instability |10| . 

4 Analysis of linear stochastic system: subcritical case 

In this section we discuss the circumstances under which the magnetic field of initially zero 
amplitude can experience sustained growth even in the subcritical case (kinematic regime). The 
linear approximation of the stochastic system (0),® near equilibrium point = 0, B^p = 
can be written as 





-e -5 




Bj. 




1 


A = 




, B = 




, a = 






-g -e 




Bip 








a{B) = aiBl + a2 = ai{B,SB) + a2, S = 
Here W is the standard Wiener process. 





1 



(22) 
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The important statistical characteristics of this system are the second moments, since they 
represent the energy of the magnetic field (see, for example jl()(in]). The second moments matrix 
M = ii^(BB''') is governed by system 

^ = AM + MA"^ + 2(o-itr(MS) + cr2)aa'^. (23) 

There exists a variety of formal derivations of the equations l\2'6^ . One of them can be found in 
Appendix A. It is also convenient to represent the matrix M in the form 



M 



mi 1712 
m2 ms 



(24) 



where mi{t) = E{B'^), m2{t) = E{Bj.B^), m^{t) = E{B'^), and E{-) denotes the expectation. 
4.1 Additive noise 

In this section we study how the additive noise can amplify the magnetic energy 

EB{t) ^ m3(i) = E{Bl) (25) 

during the kinematic stage. Without multiplicative noise (ai = 0) the system ()2U() takes the 
form 

§^AB + V2r,.^. (26) 
at at 

It follows from ()23() that in the subcritical case the second moments matrix M(t) converges to 
a stationary value as t ^ oo such that 

(2e2 - gS)a2 g(T2 g'^(T2 

= 2e{e^ - 95) ' "^^ = 2(^5^' = 2e{e^ - gS)' ^'^^ 

This result implies that additive noise effects, which occur naturally in the turbulent conducting 
fluids, together with differential rotation provide mechanisms which generate the large scale 
magnetic field in the subcritical case. Let us consider the dependence of the stationary value 
rris = E{B'^) on the differential rotation parameter g. When g varies from to Qcrit — ^ 

/5 the 

function m^{g) increases such that m-i{g) — > cxd as 5 — > gent = e^/(5. Note that the sensitivity 
of normal dynamical systems to additive noise near the bifurcation point is a well-known result 
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Let us discuss now the differences between a non-normal system like (|2Hj) and the normal 
one in the presence of additive noise. To understand this difference consider along with the 
non-normal system (|2(ij) the corresponding scalar normal stochastic equation 

^ = X,B + V2^2^, Xi = -e + VgS. (28) 

In the subcritical case (Ai < 0), the solution of H28() is nothing else but the classical Ornstein- 
Uhlenbeck process j3Ul I31j . It is easy to find the stationary value of its second moment m = 
E{B'^) : 

m = ^^. (29) 

To assess the significance of the effect of non-normality with respect to the sensitivity to additive 
noise, let us introduce the new parameter k = ^ that can be interpreted as a stochastic non- 
normality coefficient. If we consider the neighborhood of the bifurcation point Ai = 0, it gives 
us the measure of sensitivity of the non-normal system to noise compared to the normal one. 
By using and (^H), one can find 

k = , (30) 
2e(e + V5?) ^ ^ 

At the bifurcation point, g = Qcnt the value of this parameter is 

k = %t- (31) 

Since Qcrit ~ 1 it follows that k ~ >> 1, which shows how sensitive the non-normal 
system is compared to the equivalent normal system We note that the level of second 
moments maintained in stochastic non-normal dynamical systems associated with linearly stable 
shear flows has been discussed in 

4.2 Multiplicative noise: stochastic instability 

In this section we discuss how the average magnetic energy Es^t) ~ = E^B"^) is amplified 
by the random fluctuations of the a— parameter in the subcritical case (kinematic regime). In 
particular we find the critical value of the multiplicative noise intensity acr such that for all 
values of ai > acr the energy Esit) grows as exp(At). during the early kinematic stage. 
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Consider the system (|2()|) with multiphcative noise only 

^ = AB + y'2ai(B,SB)a^. (32) 

It is well known that the second moments matrix of the system (j2(Jj) converges to a sta- 
tionary value as t — > cxD if and only if the equilibrium B = of system (|32|) is exponen- 
tially stable in the mean square sense (EMS-stable). EMS-stability means that ^||B(t)||2 < 
i^exp(-/t)S||B(0)||, K,l>0. 

The equilibrium B = of system (|32p is EMS-stable if and only if 

a) the equilibrium B = of the deterministic system is asymptotically stable; 

b) tr(MS) < 1, where M is the second moments stationary matrix for the system 

— = AB + V2^a— . 33 

at at 

One can find the proof of this result in Appendix B. This criterion (see |38| ') reduces the linear 
stability analysis of a system with multiplicative noise to that of the second moments stationary 
matrix of the corresponding system with the additive noise only. 

For the subcritical case, it follows from this theorem that the system ()2U() is EMS-stable if 

g^a,<2e{s^-gS). (34) 

The critical value of the multiplicative noise intensity cJi is 

(Tcr = 5 . (35) 

g^ 

If 0"! > (Tcr then for any additive noise intensity a2 > the second moments of the system (|32|) 
tend to infinity as t ^ oo. This means that the random trajectories of the nonlinear system 
O,® leave the basin of attraction of zero equilibrium Br = 0, = . The understanding 
of the stabilization of this growth should involve the numerical solution of the full non-linear 
problem given by Q,®- 

It should be noted that the idea of a magnetic generation by differential rotation (the 
w— effect) is widely accepted, but the a— effect is still considered as controversial. It follows 
from the above analysis that the amplification of magnetic energy might happen even in the 
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case when the mean value of a is zero. It should be noted that this result is different from analo- 
gous in since we consider here the second moments of the random magnetic field. Certainly, 
the predictions obtained here for the simplified stochastic model should be ultimately extended 
to more realistic systems of partial differential equations ((21). However, for now we use the no— z 
model as a simplified tool to observe the noise effects on the generation of magnetic energy in 
the subcritical case. 

5 Numerical analysis: nonlinear case 

So far we have concentrated on the linear stochastic instability where the second moments grow 
exponentially without limit (kinematic regime). Obviously, in the non-linear case, if we take 
into account the backreaction which suppresses the effective dissipation j3, and the a— effect, one 
can expect an entirely different global behavior. 

In this section, we perform simulations of random trajectories of the nonlinear dynamical 
system 0,® for e = 0.1, S = 0.01, ka = = 1. In this case the critical value of the differential 
rotation parameter gcrit is 1. Since we are interested in the subcritical case, 5 < 1, we choose g 
to be 0.99. It follows from ()35|) that for this set of parameters the critical value of multiplicative 
noise intensity Gct = 2 • 10~^. Let us emphasize again that the main reason why is so small 
is because of the high level of non-normality of the dynamical system ((TJ,®. 

Fig. 3 demonstrates qualitative changes in the shapes of the probability density function 
(pdf) of B^p for fixed t as the intensity of multiplicative noise ai increases {(J2 = 0.5 • 10^"^). 
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-2 -1.5 -1 -0.5 0.5 1 1.5 2 

B 

♦ 

Fig. 3. Probability density functions for different values of the intensity of multiplicative noise. 

Here asterisks mark the positions of the equilibrium points for the deterministic system. 
One can see that there are three qualitatively different regimes depending on the value ai. For 
very small values of fJi the pdf is concentrated around the equilibrium point B^p = 0. Stochastic 
trajectories of ((TJ,© calculated by direct numerical simulations are shown in Figs. 4 and 5. 

E=0.1, 8=0.01, g=0.99, cr =2 10"^, o =0.5-10"' 

2 

1.5 
1 

0.5 


-0.5 



-0.12 -0.1 -0.08 -0.06 -0.04 -0.02 0.02 0.04 0.06 0.08 

B 

r 

Fig. 4. Stochastic trajectory escapes from zero equilibrium point and concentrates near 
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non-zero equilibrium. 




-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 

B 

r 

Fig. 5. Stochastic trajectory moves around all stationary points. 

One can see that as ai increases, the trajectories (see Fig. 4) escape from the domain of 
attraction of the zero equilibrium point and concentrate in the vicinity of the non-zero sta- 
ble equilibrium point. Further increase of cti leads to quite complicated dynamics when the 
stochastic trajectory (see Fig. 5) moves around all deterministic stationary points. 



6 Discussion and conclusions 

We have studied the stochastic amplification of large scale magnetic fields in a differentially 
rotating system in a subcritical regime and discussed the possible implications of this for the 
magnetic field in galaxies. The main purpose was to address the stochastic generation that 
cannot be explained by traditional linear eigenvalue analysis. We have chosen the simplified 
stochastic afi— dynamo model for galaxies in a thin-disk approximation and thereby concen- 
trated on the influence of additive and multiplicative noises along with non-normality on the 
amplification of the magnetic field in the subcritical case (when the dynamo number is less than 
critical). In the linear case, we have derived the equations for second moments describing the 
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magnetic energy and demonstrated the important differences between the non-normal system 
and the normal one under the influence of additive noise. For the multiplicative noise, the crite- 
ria for the stochastic instability during the early kinematic stage was established in terms of the 
critical value of noise intensity due to a-fluctuations. In the non-linear case, we have performed 
numerical simulations of non-linear stochastic differential equations for the aO— dynamo and 
found a series of noise induced phase transitions: qualitative changes in the behavior of the 
trajectories due to the increase in the noise intensity parameter. 

It should be noted that the equations as applied to galaxies are a theoretical simplifica- 

tion, but the do provide a useful framework for understanding the effect of random fluctuations. 
Our finding for the stochastic parametric instability can be straightforwardly applied to partial 
differential equations @ in which case one has to consider the coupled system for eigenmodes 
\12\ . In this respect the two equations (O,® can be regarded as a dynamical system for first 
eigenmodes (order parameters), where the influence of other degrees of freedom is parameter- 
ized by additive and multiplicative noises. We believe that our theory can also be extended and 
applied (after some modifications) to the solar dynamo in which a toroidal field is generated by 
the action of a shear flow (typical non- normal effect). 
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Appendix A 

Formal derivation of the equation for the second moments matrix 

Let aB = B(t + At) - B(t). It follows formally from ^ that aB = ABAt + vaW, where 
B = B(i), V = y/2a{B{t)) ■ a or 

B(t + At) = B + ABAt + vAT^ (A-1) 

Equation ()A-1|) implies 

Bit + Ai)B'^(t + At) = (B + ABAt + vAl^)(B^ + B"^A"^At + aW) = 

= BB^ + (ABB"^ + BB'^A'^)At + vv'^(ATy)2 + (vB"^ + Bv'^)aT4^ + 

+ (ABv^ + vB^A^)AtAT4^ + ABB^A^(At)2 (A-2) 

Taking into account the standard relations EaW = 0, E{aW)'^ = At for a Wiener process and 
£;(vv^) = 2{aitr{MS) + a2)aisJ we can find from for M(t) = E{B{t)B'^ (t)) the following 

equation 

AM = M{t + At) - M{t) = (AM + MA^ + 2{aitr{MS) + a2)aa^)At + AMA^ (At)"^ (A-3) 
Now system ()23p follows from ()A-3p immediately. 

Appendix B 

Necessity. From EMS-stability of H32() it follows that there exists a stationary second moments 
matrix M of system H20p satisfying the equation 

AM + MA"^ + 2(f7itr(MS) + o-2)aa'^ = 0. (B-1) 

One can find for ai > that the matrix 

M = M/ (^tr(MS) + (B-2) 

is a solution of the equation 

AM + MA'^ + 2(71 aa^ = 0. (B-3) 



20 



Note that M is a stationary second moments matrix for the system (jHSJ • The obvious inequaUty 

(B-4) 



,,(MS) = -4^^^ < 1 



£2 
0"1 



tr(MS) + 
proves necessity. 

Sufficiency. Consider a stationary second moments matrix M of system ()c{3|) satisfying 
equation HB-3|) and the inequahty tr(MS) < 1. It means for sufficiently smah e > a solution 
Mg of the equation 



AMe + M^A^ + 2(7iaa^ + 2el = 0, I 



1 
1 



(B-5) 



satisfies the inequality 



too. From (1^31). (TB^ it follows that 



tr(Me5) < 1 



(B-6) 



AM,, + M^A^ + 2{aitv{MeS) + a2)aLsr + 2el = 



(B-7) 



where (T2 = (Ti(1 — tr(M£S)) > 0. Formula ()B-7|1 shows that is a stationary second moments 
matrix for the stochastic system 



dB 
'dt 



dW 



d^ 



AB + y^2(ai(B, SB) + f^2)a— + 



(B-8) 



where E, is the independent two-dimensional standard Wiener process. The existence of a station- 
ary second moments matrix for the stochastic system (jB-Sp with two-dimensional nondegenerate 
additive noise of intensity e > proves EMS-stability of the system ((3H). The details of the 
proof can be found in [HH] . 
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